Additivity Principle in High-dimensional Deterministic Systems 
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The additivity principle (AP), conjectured by Bodineau and Derrida [Phys. Rev. Lett. 92, 180601 
(2004)], is discussed for the case of heat conduction in three-dimensional disordered harmonic lattices 
to consider the effects of deterministic dynamics, higher dimensionality, and different transport 
regimes, i.e., ballistic, diffusive, and anomalous transport. The cumulant generating function (CGF) 
for heat transfer is accurately calculated, and compared with the one given by the AP. In the 
diffusive regime, we find a clear agreement with the conjecture even if the system is high-dimensional. 
Surprisingly even in the anomalous regime the CGF is also well fitted by the AP. Lower dimensional 
systems are also studied and the importance of three-dimensionality for the validity is stressed. 
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Understanding the general features of currents and 
their fluctuations in transport phenomena is one of the 
main goals in nonequilibrium statistical physics. Heat 
conduction is a typical transport phenomenon, where 
one considers heat transferred through a system from a 
bath at temperature to a bath at temperature Tr. 
Transport can be classified by the system-size depen- 
dent thermal conductivity defined in the linear response 
regime. Considering a slab of width W and length iV, 
for a small AT — Tl — Tr applied across its length and 
with T ={Tl+ Tr) /2 wc define 



^ ' AT/N 



(1) 



where J is the energy current. Fourier's law implies a — 
0, while a = 1 is ballistic transport. Many dynamical 
systems with momentum conservation show anomalous 
transport (0 < a < 

In this letter, we consider the properties of current 
fluctuations beyond the linear response regime in differ- 
ent transport regimes characterized by the parameter a. 
One of the universal properties of current fluctuations 
which is expected to be valid irrespective of a, is the fluc- 
tuation theorem 0-01 ■ quantitatively connects distri- 
bution of positive and negative heat transfer and is valid 
in the far from equilibrium regime. For diffusive systems, 
where Fourier's law is satisfied (a = 0), some irnportant 
progress has been made. Bodineau and Derrida [8| made 
a remarkable conjecture, namely the additivity principle 
(AP), which enables one to compute all higher orders 
of current cumulants given just the temperature depen- 
dent thermal conductivity k{T) of a system. Bertini and 
co-workers Q introduced a macroscopic fluctuation the- 
ory (MFT) that describes the asymptotic probability of 
observing a given time-dependent local current and tem- 
perature proflle. The MFT is expected to be valid for 
a wide class of stochastic models, and the AP can be 
derived from it under the condition that the dominant 
trajectories are time- independent. From this, the suffi- 



cient condition to get AP from the viewpoint of the MFT 
is given by @] 



k{T)[k{T)TJ' < [k{T)]'[k{T)T'^ 



(2) 



However, we do not still have the necessary and suf- 
ficient condition for the AP. One of the strategies for 
finding the condition is to test the AP in different con- 
crete models. It was confirmed that the AP is consis- 
tent with the exact expressions of several orders of cur- 
rent cumulants in the symmetric simple exclusion pro- 
cess 3. Recently, the AP was numerically verified in 
another stochastic system, namely for heat transport in 
the Kipnis-Marchioro-Presutti model [ll]|, by measuring 
rare events with a sophisticated algorithm 12[. However, 
studies so far are concerned only to stochastic processes 
where not only the reservoir but the system dynamics 
also is probabilistic. We do not still understand the 
effects of deterministic system dynamics, higher dimen- 
sionality, and non-diffusive transport (a ^ 0). Hence, it 
is of general interest to consider the AP for systems with 
bulk Hamiltonian dynamics attached to stochastic ther- 
mal reservoirs. In this letter, we for the first time address 
these effects using the three-dimensional mass disordered 
harmonic crystal which was recently shown to show dif- 
ferent regimes of transport, such as ballistic {a ~ 1), 
diffusive (a = 0), and anomalous transport (0 < a < 1) 

p. 

Let Q be the heat transferred from the left reservoir 
to the system during the measurement time r, and let 
P/v((7, Tl, Tr) (g = Q/t) be the distribution of Q for the 
system with the size N . In general, the distribution at 
large r has the large deviation form TV ^ QThN{q,TL.Tii) ^ 
where hi^ is the large deviation function (LDF). The AP 
states that the LDF is given by the sum of the LDFs 
of subsystems of length n and N — n: hfq[q,TL,Tii) — 
maxT [hn{q,TL,T) -\- hjq-n{<l,T,Tji)\ . This conjecture is 
applied iteratively to break the system into a number of 
smaller pieces. Using the assumption of local equilibrium 
for the small pieces one can obtain an explicit expression 



2 



for the LDF of the full system [8|. The cumulant gener- 
ating function (CGF) defined as 



1 



^(A) = lim -log(e^^) 



(3) 



is connected to the LDF through the Legendre transfor- 
mation /i(A) = maxg [Xq + /iAr(g, Tl, Tr]J and from this 



one also gets an expression for the CGF [8|. Verifying the 
AP prediction through direct simulations of heat conduc- 
tion in a 3D crystal is extremely difficult. In this letter 
we use some recent exact results on the CGF of a har- 
monic crystal to test the AP. 

Model and Methods. — We consider a 3D cubic har- 
monic crystal with a scalar displacement field Xn on each 
lattice site n ~ 7121^3) where ni — 1,2, . . . , N and 
'^2j "-3 = 1; ■ • • W. The Hamiltonian is given by 



n 



+ y4 + ^ 2^(a^r^ - (4) 



where e denotes unit vectors in the three dimensions. We 
have set the spring constant between sites to one. Masses 
are randomly distributed as 



1 - A or 1 



(5) 



with equal probability. Two faces of the crystal, namely 
those at rii = 1 and ni — N, have fixed boundary condi- 
tions and are coupled to white noise Langevin type heat 
baths at temperatures and Tr (< Tl), respectively. In 
the transverse directions, periodic boundary condition is 
imposed. Let £ and r be the sites of left and right faces, 
namely £ = (I,n2,n3), and r ~ (N ,nl^,nl.^. The equa- 
tions of motion of the particles are then given by 



-koXn - y^(a 

e 

Sn,£h±£ -yi) ~y^, Sn,r{lXr - Vr) ■ (6) 



The noise terms at different sites are uncorrelated, while 
at a given site the noise strength is specified by the corre- 
lations {r]i{t)r]t{t')) = 2-/TLS(t-t') and {r]r{t)r]r{t')) = 
2'yTfi6{t — t') where we have set the Boltzmann constant 
to the value one. 

We assume that the initial state at i ~ is chosen 
from the steady state distribution . The heat Q flowing 
from the left reservoir into the system between the times 
t = to t = r is given hy Q ^J^e lo ^-^ ("7*^ + 
The average current J = (Q) /r in the harmonic crystal 
is given by a landauer-like formula [lo| and this gives 
the following expression for the size-dependent thermal 
conductivity defined in Eq. ([1]): 



N 



J/W'^ _ 
AT/N ~ 2ttW^ 



du;Tr[T{uj)], (7) 



A[G+TrG-Tl\u' , 
[-Mc^2 ^ -1 ^ 
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FIG. 1: The size-dependent thermal conductivity k for 
(A, fco) = (a): (0, 0), (b): (0.22, 9.0), and (c): (0.82, 0). For suf- 
ficiently large W we see respectively (a): ballistic {a = 1), (b): 
diffusive (a = 0), and (c): anomalous transport(a ~ 0.23) . 
Each data-point is obtained from simulations of Eq. ((Sjl , while 
the black dotted line is given by Eq. (0 for largest W . 



where T{uj) is the transmission matrix which describes 
transmission of phonons emitted from a site on one face 
attached to a reservoir to a site on another face, and 
is a VF^ X W'^ matrix. The Green's function G^iuj) 
is a NW"^ X NW"^ matrix given by the mass and force 
constant matrices M and K , and the self-energy matri- 
ces whose matrix elements are \^L/R\n,n'{^) — 
«7'^'5n,n"5„,£/r • The matrix Tl/r = Ivcl{T, l / r{uj)} . 

Recently it was shown 14] that, not just the current, 
but an exact expression for the full CGF of the harmonic 
crystal can be obtained in terms of T{uj) and is given by: 



IJ-HcW 



duj 
2^ 



Trlog 



T{uj)TlTrX{X + Ai3) , (8) 



where (3l = TZ\(3r = T^^ and Ap ^ (3r- (3l- One 
can verify that the current in Eq. ([7]) is given by J = 
^^ol^ \\=o- Also Eq. ^ satisfies the fluctuation theorem 
symmetry relation: /i(A) = IJ.{—X — A/3) 

We next discuss the prediction for CGF from the AP, 
which we will denote by /z^p(A). In general the CGF can 
be expressed completely as a temperature integral over 
the range [Tr,Ti] involving the single parameter k{T). 
For the harmonic case, k is independent of temperature, 
and explicit expressions can be obtained for /i^p(A) [ll| 
in terms of a single parameter k. These expressions are 
somewhat lengthy to state and hence we give them in the 
supplementary material We note that ^apW also 
satisfies the symmetry: /i(A) = A — A/3). 

The main aim of this letter is to compare the AP pre- 
diction for fj,Ap{X) with the numerical result for fiHcW 
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FIG. 2: Temperature profiles for each case in Fig.[T] Temper- 
ature is defined as Ti = I]„.„j=,(m,„i^)/W^. 



from Eq. We note that for a 3-D disordered crys- 
tal the heat current depends on the particular realiza- 
tion of disorder, however for large N and W there is 
self-averaging and sample-to-sample fluctuations become 
very small. Hence for a fixed disorder strength we get 
a unique current and n from Eq. ([T]). This value of k is 
then used to get ^apW and compared with ^hcW- 

Average heat conduction for different regimes. — In 3D 
disordered systems without pinning potentials (fco = 0), 
low-frequency extended modes with diverging phonon 
mean-free-paths exist and lead to anomalous transport. 
However, a pinning potential removes these modes and 
transport is then governed by the high-frequency ex- 
tended diffusive modes. Hence a 3-D disordered pinned 
crystal shows diffusive heat conduction. Based on the 
results of [13] we expect different regimes of transport 
and accordingly we chose the following three parame- 
ter sets for these regimes: (a) ordered unpinned lattice 
(A = 0, fco = 0) for ballistic transport {a ~ 1), (b) dis- 
ordered pinned lattice (A = 0.22, fcp = 9.0) for diffusive 
transport (a = 0), and (c) disordered unpinned lattice 
(A = 0.82, fco = 0) for anomalous transport (0 < a < 1). 
To demonstrate the different transport regimes we show 
the size-dependence of thermal conductivity in Fig. [l] 
and the typical temperature profiles in Fig. [2j The aver- 
age heat current was obtained either by direct nonequi- 
librium simulations of the Langevin equations Eq. ^ , or 
from Eq.([7l) using recursive Green's function techniques 
[isj to evaluate the transmission matrix T{uj). Agree- 
ment between the two methods is excellent. Each point 
in Fig. [1] is for one disorder realization and we fixed pa- 
rameters 7 = 1, Ti = 2.0, and Tr = 1.0. 

Remarks on Figs. [T] and [5] are in order. From Fig.[TIa) 
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FIG. 3: Comparison of the numerically evaluated CGF 
fJ-Hc{^) (points) for a 3-D harmonic crystal and the AP pre- 
dicted curve jj,Ap{\) for N = 128 and different widths. Sys- 
tem parameters are the same as in Fig. [l]with Tl ~ 2.0 and 
Tr = 0.25 . The black dotted line is a quadratic fit to fiHcW 
for the largest W. The range of A is (— /3_r,/3l). 



we see that k is independent of width and diverges lin- 
early with N implying ballistic transport. On the other 
hand in Fig. [IJb) andfTfc), we see that for small W, k 
decreases for increasing N. This implies the emergence of 
phonon localization since the system is quasi-one dimen- 
sional. For increasing W with fixed N, the data converges 
to a constant value, which implies the self-averaging ef- 
fect in disordered systems. Hence, one can get precise 3-D 
behavior for sufficiently large W. Fig. [DJb) shows diffu- 
sive transport (a = 0) for sufficiently large W and N, 
while Fig. [Ijc) shows anomalous behavior with system- 
atic power law divergence (a « 0.23). The temperature 
profile in Fig. [IJb) shows clear linear profile consistent 
with the Fourier's law. Interestingly, even the anoma- 
lous case in Fig. [D^c) shows a linear profile which is very 
different from other nonlinear systems with anomalous 
transport, such as the Fermi-Pasta-Ulam (FPU) chain 
which have nonlinear temperature profiles even for small 
temperature difference between reservoirs 0, Q • 

The CGF. — We now present results comparing 
l^Ap{^) with fiHcW ill the three different regimes . We 
obtained ^^p(A) by using k from Eq. ([7]) while fiHcW 
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was computed from Eq. ([8]). In these computations, we 
set the parameters 7 = 1.0, = 2.0, Tr = 0.25 and size 
N — 128. The results for ballistic, diffusive, and anoma- 
lous cases are shown in Figs. [2Ia), EJb), andEJc) respec- 
tively. Note that both hhcW in Eq.Q and ^apW in 
[l6j exactly satisfy /x(0) = /z(A/3) — 0. For ballistic case 
(a), we see deviations from AP curve irrespective of W 
as expected. On the other hand, in Fig.[3l[b), we see that 
the agreement between /i^p(A) and fJ.HcW improves for 
increasing W. For small W where localization effect is 
dominant, there are clear deviation from AP curve. For 
other cases with larger temperature differences, we ob- 
tained agreement with the AP for sufficiently large sys- 
tem size [1^. We should also note that the present sit- 
uation does not satisfy the sufficient condition from the 
MFT ©, since the thermal conductivity is independent 
of temperature. Hence, this model extends the sufficiency 
condition for the AP. Surprisingly, for the case of anoma- 
lous transport, we see from Fig.[3l^c) behaviour similar to 
the diffusive case, with clear verification of AP at suffi- 
ciently large W. Note that the original theory of the AP 
conjecture @ assumes diffusive transport. 

The degree of coincidence seen in Fig. [3] is now quan- 
titatively discussed. We note that harmonic lattices, in 
not only 3D, but also in ID and 2D can show diverg- 
ing thermal conductivity. For instance, ID disordered 
harmonic chains with open boundary condition show di- 
verging thermal conductivity with the power a = 1/2 
[TtI] . Then, one may ask if low dimensional anomalous 
transport satisfies AP or not. Hence, in addition to 3D 
cases in Fig. [U we also discuss low dimensional harmonic 
systems. We define the following quantity 



where A* is the value of A which minimizes /iyip(A). As 
seen in Fig. [Sj the deviation becomes maximum at the 
minimum value of the CGF. Hence, the function S quan- 
titatively estimates the degree of discrepancy. In Fig. SI 
we show i5 as a function of W for the three cases in 3D. 
Systematic approach to the AP is seen on increasing W 
for both diffusive and anomalous cases. In the inset, ID 
and 2D results for S are shown. We consider the system 
size N for ID and N x N for 2D with open boundary 
condition, and hence, the x-axis is N (not W). For both 
ID and 2D, results for one realization of random mass 
are shown. Contrary to what happens in 31?, in low di- 
mensions we see no sign of decay of 6, and it remains 
almost constant value. This implies that the coincidence 
for anomalous transport is true only for 3D systems. 

Discussion. — We have discussed the additivity prin- 
ciple in high-dimensional deterministic systems and con- 
sidered the effects of deterministic dynamics, dimension- 
ality, and non-diffusive transport (a 7^ 0). The AP was 
originally proposed for ID diffusive systems. Our main 
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FIG. 4: The quantitative estimation of the degree of agree- 
ment of the CGF with AP. The inset shows the results for 
low-dimensional anomalous cases. {Tl,Tr) — (2.0,0.25). 



result is to show the validity of AP in a 3-D Hamiltonian 
system in both the diffusive and anomalous regimes. 

In disordered harmonic crystals in ID and 2D we find 
that AP is not satisfied for the anomalous cases where 
< a < 1. The major difference between 3D and lower 
dimensions is that in 3-D only a small fraction of the nor- 
mal modes are localized while in lower dimensions, most 
of the modes are localized 



13l |18[ . Hence it is expected 



that there is no local equilibration in low-dimensions. 
Then neither the MFT nor the AP are satisfied in these 
regimes. Our study thus suggests that only the assump- 
tions of (i) local equilibration and (ii) a small current 
(requiring a < 1) are necessary for the validity of AP. 
The realization of AP in diffusive and anomalous cases 
should be related to the time-independence in dominant 
trajectories in the MFT |9| . A verification for this would 
be an important future problem, but it is not possible 
within the present approach, and would require simula- 
tions with algorithms as in |12| . 

The mechanism of anomalous transport in disordered 
harmonic lattices has some crucial differences from that 
in low-dimensional nonlinear systems such as FPU chain 
where Levy type of energy diffusion is seen [3l and tem- 
perature profile is always nonlinear. For disordered har- 
monic crystals the linear temperature profiles suggest 
that a local response relation j(x) — —ndT/dx is always 
valid, but with a size-dependent k in the anomalous case. 
On the other hand levy type of diffusion seen in one- 
dimensional nonlinear models means that the response 
is non-local (in space) and the AP may need modifica- 
tion. Hence, it will be of great interest to understand 
the general features of CGF in low-dimensional nonlin- 
ear system. 

KS was supported by MEXT (23740289). AD thanks 
DST for support through the Swarnajayanti fellowship. 
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Supplementary Material for 
Additivity Principle in High-dimensional Deterministic Systems 



Cumulant generating Junction from the additivity 
principle. — We derive the explicit expression of cumu- 
lant generating function (CGF) from the additivity prin- 
ciple (AP) for the regime -/3ij < A < 13 l- The AP as- 
sumes that the large deviation function (LDF) for the 
system with the size N is given by [l| 



so that the optimal profile remains monotonic. Then 



dT{x) 
dx 



-Nqy/1 + AKT'^k 



(13) 



hN{q,TL,TR) 



■ min — / dx' 

T{x) N Jo 



Nq + k{T) 



dT{x) 
dx 



(10) 



where k{T) = W^k{T). Variational problem in Eq.([TO]) 
is reduced to finding the optimal profile T(x) satisfying 



dT{x) 
dx 



(iVq)2 [1 + 4KT'^k{T)] 
k^{T) 



(11) 



where the function K is determined from the boundary 
condition r(0) = Tl and T(l) = Tr. From now on, 
we consider the case of temperature- independent thermal 
conductivity k(T) = k, which is the case for harmonic 
crystals. We then derive the explicit expression of the 
CGF. The CGF is given by the Legendre transformation: 



/i(A) = max [qX 

9 



hN{q,TL,Tii)] 



(12) 



In this case, from the Legendre transformation, the CGF 
is given by [l[ 



A = 



K 


/ rIT 


k 


2 






Vi- 


^AKT^k\ 




r 


dT-^ 




1 


- 1 






h AKT'^k 



(14) 
(15) 



This expression is valid for A_ < A < A+ where A-j- ~ 
(Pl - I3r± - /3i)/2- To go beyond this regime of 
A, we consider a nonmonotonic optimal profile given by 



dTjx) 
dx 




VI- 











0<x<x^, (16) 

Xc < X < 1 , 



Suppose Tl > Tr, and the deviations are not too large 



where Xc satisfies dT{xc) / dx 



— implying 1 -|- 
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0. In these cases, the CGF is given by 



N 

T(x 



T{x,) 



dT- 



Tl 

dT- 
dT 

I 

dT 



1 



(17) 



2T2 
1 

272 



1± 



1 



-1± 



1 



(18) 



These expressions are valid for the regime A+ < A < /3i 
and — /3_R < A < A_, respectively. 

Simphfying Eqs. ([14]), (US]), ^ and (HH) is straight- 
forward. By evaluating the integrations we get the fol- 
lowing explicit expressions: 



MAP (A) 



iN 



log 



y/l+4kKTl + y/4kKTl 

•••A_ < A < A+, 



(19) 



4Af 



-/3r<A<A_ 
X+<X<I3l ' 



where 9a (a — L,R) is given by 



sm9a = Vik\K\l2- 
The function K is given by 



(20) 
(21) 



KiX,TL,TR) = — {(iL 
Idk L 



Pi 



PI -PI 

PL~PI- 



2xy 

2 



2(/3i 



5|) 



2A 



(22) 



Effects of small system sizes and large temperature 
differences. — Since the AP is derived under the assump- 
tion of local equilibrium as in Eg. dlOp . one expects that 
smaller systems show larger deviation from the AP due 
to violation of local equilibrium. To see this, we con- 
sider 3-D systems with cubic structure, and calculate 
the CGF for N = 4,16,32 with {Tl.Tr) = (2.0,0.25). 
In FiglSl the CGF and the deviation 5 defined in the 
main text are shown for these cases. As expected, = 4 
clearly shows deviation from the AP curve. As increasing 
N , the deviation 5 decreases. 

We next consider the effect of large temperature dif- 
ferences for large systems. For this we obtain the CGF 
for many temperature sets at fixed large N . We note the 
following exact scaling relation that holds for the CGF : 



(23) 



where v is an arbitrary real number. This holds for both 
^AP and ^HC- From this relation, it is clear that if 
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FIG. 5: The CGF in small systems. The system is a cubic 
structure A*'^. 



agreement between the formulas is seen for some tem- 
perature difference AT = Tl - Tr {Tl > Tr) and av- 
erage temperature T — {Tl -t- Tr)/2, then it also holds 
for i/AT and vT . Hence the correct relevant parameter 
is the ratio AT/T. In FiglHJ we show the CGF of many 
temperature sets with N = 128. Points are fJ-Hc/W^, 
while solid lines are ^ap/W^, and the regime of A is 
A G [— /3_r,/3l]. We see that for the ordered lattice the 
disagreement with AP becomes large on increasing AT/T 
while for the disordered case AP is always satisfied. FiglT] 
clearly shows that in sufficiently large systems, the AP 
is accurate over a large range of AT/T and the agree- 
ment is better for smaller AT/T. The relevance of the 
parameter AT/T can be roughly understood by consid- 
ering the criterion for local thermal equilibrium. If the 
typical mean free path of the heat carriers (phonons here) 
is denoted by £p then the condition for local equilibrium 
is ip\dT/dx\/T « 1 [2] or AT/T « L/£p. For the or- 
dered ballsitic case £p ^ L and so we require AT/T << 1. 
On the other hand for disordered systems ip is finite and 
hence for sufiicienlty large size L the condition for local 
equilibrium is always satisfied for any given AT/T. We 
also note that in diffusive systems the condition for local 
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FIG. 7: Deviations from the AP in various temperature 
regime with fixed temperature difference AT = Tl~Tr — 1.0. 



equilibrium also ensures that temperature profiles do not 
show any jumps at the boundaries. 

Cumulant generating function for ordered harmonic 
crystal. — We here give the simplified expression of 
fiHcW ordered harmonic crystal (A = 0). The sim- 



plified expression for fJ,HcW is much more computation- 
ally efficient than using the recursive Green's function 
technique for solving the CGF. 

Each lattice is labelled by the vector n = nie.i+n2e2 + 
^363 where ni — I, N and ^2, ri3 = 1, W. Let C/^"^ 
be a matrix which acts only on the components in the 
CQ-direction. We introduce the orthogonal matrix 



U 



- 1(1) 



(24) 

a = 2,3. (25) 



Using this matrix, scalar variable Xn defined on the site 
n is transformed as 



cos I — — — I cos 



27r 7713713 
W 



(26) 



where the vector m stands for m = (mi, rn^), {mi, m^ = 
I, W). By this transformation, the Hamiltonian with 
A = is transformed into independent chains of 

the form H = E„, + ^i?' -^'"^^ 



7i \ ,n -|^ 



/2, 



where „/ is the m-mode representation of the force 
matrix K: 



K 



(m) 




4[sin^(^)+sin^(^)] 
• • • ni = n'l 



(27) 



= 1 



In addition, noise terms preserve the correlations 

{v[T^{t)^^^'\t')) = 2^n5,^^,,drr^^r^,S{t~t'), (28) 
UTHt)^lT'^(*'^^ = ^lTRSr,y^Sm,m'Sit~t'). (29) 

Hence, not only Hamiltonian but also Langevin equations 
are decomposed into independent Langevin dynam- 
ics. 

We now consider the transmission matrix in terms of 
which the CGF can be written. The transmission matrix 
is diagonalizcd into the m-mode 



2,. ,21(5+ ("1)12 



1,N 



(30) 



with the Green's function given by the inverse of the 
tridiagonal matrix 



«7'^'5„j^„'^((5„i,i -I- (5„i,Ar). (31) 



The expression of is readily obtained, hence, we 

finally get the explicit formula of fJ-HcW for ordered har- 
monic crystal as follows. 



duj log 



47^^^ sin^ I 



-TLTRXiX + Al3) , (32) 



where 

2 + fco - , „ r . 2 /7I"m2\ . 2 /7I'TO3\ 
Arri = [{l - J^LO"^) COS 0m - '2ijL>j]sm{N 9m) 

+ (1 + 7^0;^) sin6im cos(iV(9^) . 
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